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Abstract 

Isotope shifts in dielectronic recombination spectra were studied for Li-like A Nd 57+ ions with 
A=142 and A=150. From the displacement of resonance positions energy shifts Ji? 2,150 (2s — 
2 Pl/2 ) = 40.2(3)(6) meV ((stat)(sys)) and <5£ 142 > 150 (2s - 2p 3/2 ) = 42.3(12)(20) meV of 2s - 2 Pj 
transitions were deduced. An evaluation of these values within a full QED treatment yields a change 
in the mean-square charge radius of 142 450^ r 2^ = -1.36(1)(3) fm 2 . The approach is conceptually 
new and combines the advantage of a simple atomic structure with high sensitivity to nuclear size. 

PACS numbers: 34.80.Lx, 31.30.Jv, 31.30.Gs, 21.10.Ft 
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The extent and shape of the proton distribution are basic properties of the atomic nu- 
cleus and reflect the interplay of forces that act between the nuclear constituents, protons 
and neutrons. The gross properties of nuclear charge radii across the chart of nuclides — 
and also the fine details within isotope or isotone chains — provide primary observables for 
testing nuclear models . Charge radii, nuclear shape and hyperfine structure studies 

with classical methods such as electron scattering and x-ray spectroscopy of electronic and 
muonic atoms are essentially restricted to stable nuclei. Thus those isotopes far off the valley 
of stability cannot be addressed where the predictive power of nuclear models is poor, and 
interesting phenomena, like neutron or proton skins and shape effects occur. This limitation 
is overcome in optical isotope shift (OIS) measurements which are well suited to study long 
isotopic chains including exotic nuclides 4|. OIS is highly sensitive and possesses an enor- 
mous experimental precision. Since measurements of OIS are usually limited to low charge 
states, the interpretation of the data is, however, hampered by atomic many-body effects. 
The difficulty of a reliable theoretical description of these complex electronic configurations 
results in large uncertainties due to specific mass shift contributions and in the calculation 
of the change of the electronic density at the site of the nucleus. 

Detailed knowledge about nuclear properties and their influence on electron binding is essen- 
tial for many applications in modern atomic physics research, in quantum chemistry [5| and 
even in the determination of chemical abundances of stellar objects [6]. This is particularly 
true for high-precision tests of fundamental interactions like atomic parity non-conservation 
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3, 0, Ei] or quantum electrodynamics (QED) in strong electromagnetic fields [lO, 
Nuclear effects play a critical role in such tests and limit their explanatory power. The con- 
straints are increasingly more severe for heavier elements since the scaling of nuclear effects 



with the nuclear charge Z exceeds the one of the sought-after contributions 



Radii compilations such as 
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, ll6j are widely used for reference. A closer look into the 
publications that form the basis of the compilations very often reveals inconsistencies. For 
instance, for the present case of neodymium, more than 20 publications from all four classi- 
cal techniques can be found (for references see (4, [l5, [if]]). In 17] a decreasing mean-square 
charge radius for a change from isotope 142 Nd to 150 Nd is reported, whereas in 18| a radius 
increase and a value of only 142 > l50 S(r 2 ) = -0.22 fm 2 are reported. While these values are 
considered as outliers the bulk of the data still covers a range from about 142 > 150 <5(r 2 ) = -1.20 
fm 2 to -1.36 fm 2 . 
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In view of the wide range of applications and the scatter of the experimental data alter- 
native methods for a reliable determination of charge radii are highly desirable. In this 
Letter we present a conceptually different experimental access to charge radii changes that 
is particularly well suited for the investigation of heavy stable and also unstable nuclides. 
The method is based on the storage ring measurement of isotope shifts (IS) 5E A ' A ' in the 
spectrum of resonant electron-ion recombination of heavy few-electron ions. The technique 
is demonstrated for the case of the two stable even-even isotopes 142 Nd 57+ and 150 Nd 57+ of 
Li-like neodymium. Storage rings provide very clean experimental conditions since isotopi- 
cally pure, quasi-monoenergetic ("cooled") beams in a single charge state are available, 
experiments were performed at the electron cooler of the heavy ion storage ring ESR 
at GSI in Darmstadt, Germany. The benefit of using few-electron ions for IS studies was 
previously shown employing classical spectroscopy at electron-beam ion traps (e.g., 0, Q)- 
Due to the simplicity of the atomic configurations, the interpretation of the data is clear 
and without ambiguity. For the 2s — > 2p transitions of Li-like ions the electronic part can 
be treated theoretically with high accuracy. Many-body and mass effects are small and can 
be reliably accounted for. 

Resonant electron-ion recombination — also termed dielectronic recombination (DR) — can be 
viewed as a two-step process. The first step of DR is called dielectronic capture (DC) and 
is time-reverse to autoionization. Free electrons that posses matching kinetic energy can 
recombine resonantly and excite a bound electron. If — in a second step — DC is succeeded 
by emission of photons, DR is complete. For the 2s — > 2pj (J = 1/2,3/2) excitations of 
A Nd 57+ DR can be described by 

e~+ A Nd 57+ (ls 2 2s 1/2 ) -> A Nd 56+ (ls 2 2 Pj nl f )** 

-> A Nd 56+ + photons. (1) 

In Nd 57+ , DC can proceed via Rydberg states with principal quantum numbers n > 18 for 
the 2s — > 2pi/2 core excitations and n > 8 for 2s — > 2p 3 / 2 , respectively (compare Figs. [1] 



and . With increasing electron-ion co 



the series limits n — > oo are formed 
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lision energy, series of Rydberg resonances up to 



221 ] . For different isotopes, the whole resonance 



series exhibits the same energy shift, i.e., the one of the corresponding core excitation. The 
captured loosely bound Rydberg electron itself is almost unaffected by the small variations 
of the nuclear potential and its contribution to the isotope shift can be safely neglected. 
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FIG. 1: Dielectronic recombination of the Li-like neodymium isotopes 142 Nd 57+ (orange line) and 
150 Nd 57+ (black line) in the energy range of the Is 2 2pi / 2 18lf resonance groups. The labels indicate 
the individual fine structure components j' of the n = 18 Rydberg electron. 

In recent DR experiments at the storage rings CRYRING in Stockholm and TSR in Heidel- 
berg, the high resolution that can be achieved for exceptionally low-energy resonances was 



exploited to explore magnetic hyperfine effects in Cu-like 207 Pb 53+ 23f| and in Li-like 45 Sc 18+ 



24j . Low-lying resonance are well known for selected ions of intermediate charge state but 
are unlikely for heavy few-electron ions. The Li-like ions under study in this Letter are ideal 
candidates for precision DR-IS experiments and provide an explicit and straightforward ac- 
cess to nuclear parameters. The 2s electronic wave function and the nucleus posses large 
mutual overlap leading to strong nuclear size contributions. In any element, DR resonances 
with excitations 2s — > 2pi/ 2 and in some cases additional resonances belonging to 2s — 2p 3 / 2 
transitions can be found in the favorable energy range below a few ten eV. 
Detailed descriptions of the recombination set-up at the electron cooler of the ESR are given 



elsewhere 
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251 ] . Here, we concentrate on the issues specific for the IS experiment with 
Li-like A Nd 57+ : Momentum selected A Nd 57+ ions with A=142 or A=150 were injected into 
the storage ring ESR and stored at an energy of 56.3 MeV/u. In order to minimize lateral 
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beam size and energy spread electron cooling was applied. For both isotopes an identical 
space-charge corrected cooling voltage of Uc = 30899 V was chosen, equivalent to an ion 
velocity /3j = 0.3328 in units of the speed of light. In addition to its normal function as 
a beam cooling device, we operated the electron cooler as a target of free electrons. After 
a series of recombination measurements with the reference isotope A=142, the operation 
parameters of the accelerator and the storage ring were scaled with the known mass ratio 
to the new isotope A=150. Special care was taken to ensure the same alignment of the ion 
beam with respect to the electron beam. Non-zero collision energies were introduced by ap- 
plying a sequence of swiftly modulated voltage steps Ud to a cylindrical drift tube located in 
the overlap region of the two beams. The relative collision energy E cm in the center-of-mass 
system can then be easily inferred from the Lorentz factors 7$ — 1 = ellc / (m e c 2 ) of the ion 
and 7 e — 1 = e(Uc + Ud) / (m 0ie c 2 ) of the electrons. Here, mo, e is the electron mass. With 
—2.5 kV <Ud< +2.5 kV used in the present experiment, collision energies of < E cm < 50 
eV were covered. Recombined ions were separated from the primary beam in the next ESR 
bending magnet and were registered by a multi-wire proportional chamber. A typical volt- 
age sequence contained measurements at ~4000 different voltages of 33 ms duration each. 
Intermediate cooling steps of the same duration guaranteed a constant ion energy and beam 
quality. A constant electron current I e i = 100 mA was applied, corresponding to a density 
n e — 3.1 ■ 10 6 cm -3 at cooling. Comparable numbers of ions of typically about — 2 ■ 10 7 
for both isotopes were available in the ring. Recombination rate coefficients a were derived 
by normalizing the count rate using the primary beam intensities jzj. 
Isotope shift data were taken for three different energy ranges, — 3.8 eV, 12 — 24 eV 
and 24 — 42 eV with nominal step widths of 1 meV, 3 meV and 4 meV, respectively. 
With these settings, resonance groups A Nd 56+ (ls 2 2pi/% nlj>)** with n = 18... 21 and 
A Nd 56+ (ls 2 2p 3 / 2 8lj/)** were covered. The shift of resonance positions for the two isotopes 
is immediately evident from the spectra in Figs. [1] and [2] and is clearly visible over the full 
energy range. For the extraction of the IS values, maxima, minima and inflection points 
of the resonance spectra were used. By evaluating the shift of these characteristic points, 
errors due to normalization and background subtraction are minimized, and an analysis is 
achieved that is independent of supplementary input on the DR process and independent 
of the experimental response function. First and second derivative spectra (Fig. [3]) as well 
as smoothed versions of the original spectrum were obtained from local fits of third-order 
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FIG. 2: Dielectronic recombination of A=142 (orange) and A=150 (black) Li-like A Nd 57+ . The 
labels n indicate the principal quantum number of the Rydberg electron of the Is 2 2p]/ 2 nly config- 
urations and the bold labels j' denote the fine structure components of the Is 2 2p 3/ / 2 8lj' resonances. 

polynomials for every point i of the spectrum taking into account m adjacent values on 
both sides of %. Full minimizations are performed accounting for a non-equidistant energy 
axis and for the error bars in both coordinates. The procdure is similar to the Savitzky 
and Golay smoothing algorithm 261 ] that retains fine details in the spectrum much better 



than averaging or re-binning 26j|. For the analysis m=35 was used and no significant bias 
for m=25 or m=45 was found. For the IS determination only those roots of the derivatives 
were included that are significantly above the noise level. For the lowest energy range 7 
independent IS data sets were obtained, 3 for the middle one and 1 for the energy range 
24 - 42 eV yielding in total 154 characteristic values for the 2s — 2pi/ 2 IS and 45 for the 
2s - 2p 3/2 IS. The average ISs are 5E 142 > 150 (2s - 2p 1/2 ) = 40.2(3) (6) meV and a slightly 
higher value of 5E lA2,150 (2s — 2p 3 / 2 ) — 42.3(12)(20) meV. Here, first and second parentheses 
denote statistical and systematic errors. Although statistically not significant, the experi- 
mentally observed small difference in the ISs is expected, since in the relativistic case the 
Pi/2-electron has a finite overlap with the nucleus. The main sources of systematic errors 
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FIG. 3: First (upper panel) and second (lower panel) derivatives of the DR resonance spectra of 
the A Nd 56+ ls 2 2p]/2 181 ji group. The energy range is the same as in Fig. [TJ 

stem from "slope" uncertainties, i.e., due to still remaining imperfections in background 
subtraction and normalization and from the finite digitalization of the energy steps. As a 
consequence of the merged-beams geometry, only minor contributions of below 0.2 meV for 
the 2pi/2 value and of ~ 0.6 meV for the 2p 3 / 2 IS arise from a potential beam misalignment 
or voltage calibration uncertainties of Uq and U4. 

The IS values were evaluated within two different fully relativistic atomic structure methods, 
namely, the multiconfiguration Dirac-Fock (MCDF) and the configuration interaction Dirac- 



Fock-Sturmian methods (see [2l| and references therein). The nucleus was modeled with a 
two-parameter Fermi distribution p(r) = iV(l + exp[(r — c)/a]) _1 , where c is the half-density 
radius, N is a normalization constant, and t= 4a In 3 the skin thickness Q. IS measure- 
ments are basically insensitive to details of the charge distribution as long as the same 



root-mean-square (rms) radius is reproduced (e.g. 
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271] ). We employed t = 2.3 fm for 
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both isotopes in accordance to the established value used for heavy ions iJ, [15|, l20l. |27] . 
The semi-magic isotope 142 Nd (82 neutrons) served as a reference using the rms radius of 
(r 2 ) 1 / 2 = 4.9118 fm from the most recent compilation of Angeli 16] . The radii values of An- 
geli were obtained from a combined analysis of comprehensive experimental input data from 
the four classical methods. The neutron-rich isotope 150 Nd (N=90) is strongly deformed 
yielding a pronounced increase of the charge radius beyond the liquid drop model 4^. 
For the 2s — 2p transitions of A Nd 57+ the total mass shift A=142— A=150 is as small as 
~4 % of the IS and sums up to -1.63 meV for the 2s — > 2p x / 2 excitation and to -1.80 meV 
for 2s — > 2^3/2, respectively. These values for the mass shift were obtained taking into 
account relativistic and QED recoil contributions 2l|, |28(: For the 2s — 2p 1 / 2 (2s — 2p 3 / 2 ) 



transition the mass shift comprises of -2.44 meV (-2.53 meV) from averaging of the non- 
relativistic recoil operator with relativistic multiconfiguration wave functions, 1.14 meV (1.03 
meV) from the relativistic recoil operator, and -0.33 meV (-0.30 meV) from the QED recoil 
effect. For the determination of S(r 2 ), influences of nuclear size (NS) variations on Coulomb 
electron correlation, on the Breit interaction as well as on QED contributions were taken 
into account. Similar to OIS Q|, to a good approximation the IS depends linearly on S(r 2 ), 
yielding a slope of -31.1 meV/fm 2 for the 2pi/ 2 and -32.3 meV/fm 2 for the 2p 3 / 2 transition, 
respectively. A further small correction that is often neglected |4J arises from nuclear po- 
larization (NP). It is sizeable for deformed nuclei with high nuclear transition probabilities 
B(tL) such as 150 Nd. Using the analytic formulae of 29| we estimate for the dominant 
contribution from the first excited 2+ state at 130.21 keV (B(E2f)=2.760 e 2 b 2 , [30]) a value 
of -0.3 meV. Finally, from the NP corrected value of 5E U2 > 150 (2s - 2p l/2 ) = 40.5(3) (6) meV, 
a change of 142 > 150 <5( r 2 } = -1.36(1) (3) fm 2 is derived. Consistently, we obtain 142 ' 150 <5<y 2 } = 
-1.38(4) (7) fm 2 from the 2p 3 / 2 resonances but with significantly larger error bars. The sys- 
tematic errors include additional theoretical uncertainties of 0.01 fm 2 , predominantly from 
NP. We have added the full NP contribution to the error balance, since NP still lacks ex- 
perimental verification. Our result of -1.36(1) (3) fm 2 is somewhat larger than the average 
experimental value of 142 ' 150 <5(r 2 ) = -1.291(6) fm 2 from the combined analysis in Angeli 16]. 
However, our values obtained with this alternative method are in good agreement with the 
upper end of the individual data that build the basis of this evaluation (see discussion on 
the first page) and support a larger increase of the radius from A=142 to A=150. 
In this Letter a novel method for nuclear size studies was developed and its general feasi- 
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bility was demonstrated for heavy three-electron isotopes. These simple systems allow for 
a reliable analysis of the IS data within state-of-the-art structure calculations taking into 
account relativistic and QED contributions. ISs could be measured with <1 meV accuracy. 
For heavier ions a similar experimental precision can be retained. Due to the Z 5 to Z 6 scaling 
of nuclear size effects, for the heaviest ions ISs are of the order of 100 meV [31( and hence 
the sensitivity of DR to the charge distribution is significantly enhanced. In contrast to OIS, 
isoelectronic studies of different elements are possible in order to disentangle atomic and nu- 
clear contributions. High experimental resolution and large resonant atomic cross sections 
enable future DR-IS experiments with unstable isotopes or long-lived nuclear isomers. At 
GSI, such exotic highly charged ions can be synthesized in the fragment separator FRS 32] 
and subsequently investigated in the ESR. Thus, long chains with constant neutron number 
(isotone shifts) become accessible as well. At the upcoming Facility for Antiproton and Ion 
Research (FAIR) [33j high production yields of radioactive ions and a dedicated electron 
target at the New Experimental Storage Ring (NESR), will furthermore increase accuracy, 
versatility and expand the scope of application. 
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